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INTRODUCTION 


Pulse  propagation  and  wave  propagation  comprise  an  important  class 
of  problems  in  the  study  of  the  dynamics  of  large  lattice  structures. 

The  study  of  pulse  and  wave  propagation  has  applications  in  dynamic 
failure,  control,  and  nondestructive  evaluation. 

In  this  investigation,  nondispersive  pulse  propagation  in  a  simple 
one-dimensional  lattice  structure  is  analyzed,  using  both  the  pulse 
summation  method  and  the  wave-mode  coordinate  method.  It  is  shown  that 
the  pulse  summation  method  (a  time  domain  method)  and  the  wave-mode 
coordinate  method  (a  frequency  domain  method)  give  identical  results, 
and  that  both  methods  account  for  the  existence  of  equivalent  paths  in 
the  lattice  structure.  Some  examples  of  equivalent  paths  are  given.  Also 
some  recommendations  are  made  for  possible  extensions  of  the  analysis 
presented  here. 


ANALYSIS 


Lattice  Definition  and  Problem  Statement 


A  one-dimensional  lattice  structure  consisting  of  two  segments  and 
one  joint  is  shown  in  Fig.  1.  It  is  assumed  that  the  joint  is  rigid 
and  massless,  and  that  the  extent  of  the  joint  in  the  x-direction  is 
small  in  comparison  with  or  Segment  one  has  elastic  modulus  E^, 

cross-sectional  area  A^,  and  mass  density  p^.  Segment  two  has  elastic 
modulus  E2,  cross-sectional  area  A2,  and  mass  density  It  is  assumed 

that  longitudinal  wave  propagation  in  each  segment  is  governed  by  the 
classical  one-dimensional  wave  equation.  Therefore,  disturbances  in 
longitudinal  force  or  longitudinal  displacement  propagate  nondispersively 
in  each  segment,  and  a  longitudinal  force  or  displacement  pulse  introduced 
into  either  segment  will  maintain  its  shape  as  it  propagates.  The  velocity 
of  pulse  propagation  in  segment  one  is 


and  the  velocity  of  pulse  propagation  in  segment  two  is 


°2  “V  P 


The  characteristic  transit  time  required  for  a  pulse  to  traverse  the 


length  of  segment  one  is 


T1  ~  c, 


m 


and  the  characteristic  transit  time  required  for  a  pulse  to  traverse  the 


length  of  segment  two  is 


(4) 


The  lengths  and  Z^  are  defined  in  Fig.  1.  It  is  assumed  here  that  the 
characteristic  transit  time  for  segment  one  is  equal  to  the  characteristic 
transit  time  for  segment  two,  or 


=  =  T 


(5) 


The  problem  considered  here  is  the  following, 
is  applied  to  the  joint  as  shown  in  Fig.  1. 


force 


is  an  impulse  of  the  form 


A  longitudinal  force 
It  is  assumed  that  the 


(6) 


It  is  desired  to  find  the  resulting  longitudinal  force  F^(t)  at  point  1, 
the  left-hand  end  of  segment  one.  Note  that  if  point  1  is  a  free  end,  then 
F^(t)  =0,  directly  from  the  boundary  condition  at  a  free  end. 


Pulse  Summation  Solution 

Using  the  pulse  summation  method  described  in  [1],  the  following 
solution  may  be  obtained  for  F^(t): 


min(K  ,K?) 

E  11  N1(K1,K2,L1)S1(K1,K2,L1) 

L!=l 


n  min(K1 ,K_-1) 

+  E  E  N  (K  ,K  ,L  )S  (K  ,K  ,L  ) 


6(t  -  (2n  +  1)T) 


where 


The  quantities  and  are  the  (displacement)  reflection  coefficients 


at  the  joint,  and  the  quantities  t^  ant  t^  are  the  (displacement)  transmission 


coefficients  at  the  joint.  The  coefficient  r^  is  the  reflection  coefficient 


for  a  pulse  which  arrives  at  the  joint  from  segment  1,  and  the  coefficient 


r 2  is  the  reflection  coefficient  for  a  pulse  which  arrives  at  the  joint 


from  segment  2.  The  coefficient  t^  is  the  transmission  coefficient  for  a 


pulse  entering  segment  1  from  segment  2,  and  the  coefficient  t^  is  the 


transmission  coefficient  for  a  pulse  entering  segment  2  from  segment  1. 


The  quantities  r^  and  r^  are  the  (displacement)  reflection  coefficients  at 


the  left-hand  boundary  of  the  lattice  (point  1)  and  the  right-hand  boundary 


of  the  lattice  (point  4),  respectively.  The  reflection  coefficient  r 


may  be  determined  from  the  boundary  conditions  at  point  1.  If,  for  example, 


point  1  is  a  fixed  end,  r^  =  -1,  and  if  point  1  is  a  free  end,  r^  =  1. 


Similarly,  the  reflection  coefficient  r^  may  be  determined  from  the  boundary 


conditions  at  point  4. 


The  methods  used  in  the  derivation  of  eqn.  (7)  are  discussed  in  detail 


in  [1].  Eqn.  (7)  is  a  slightly  corrected  form  of  eqn.  (A90)  in  [1],  and  it 


corresponds  to  the  sum  of  cases  I  and  III  defined  in  Appendix  A  of  [1] 


Eqn.  (7)  consists  of  an  infinite  series  of  impulses  which  are  delayed  by 


odd  multiples  of  the  characteristic  transit  time  t.  The  quantities  S,  S 


and  S2,  which  contribute  to  the  amplitudes  of  the  impulses,  consist  of 


powers  of  the  reflection  and  transmission  coefficients.  The  quantities 


and  ^  are  numerical  coefficients  which  will  be  discussed  and  interpreted 


subsequently.  Writing  out  the  first  few  terms  of  eqn.  (7)  gives 


^  R 

Fi(t)  =  (r,  +  rJ  (1  '  ro) 


6(t  -  T )  [  1  ]  +  S(t  -  31)1^^  +  t1r3] 
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»  f 


<  t  r 


r  2  2 


+  5(t  -  5t>|Vl  +  VlVs  +  VlVs  +  Clr2r|] 


+  $(t  -  It) 


r  1 3  .  „ 

l_r0ri  +  Vi 


t2r2r3  +  2rOrltlt2r3 


2  2  2  2  2  2  3 
+  rOrlClr3  +  r0ritlr2r3  +  rOtlt2r3  +  tlr2r3 


+  6 (t  -  9t) 


"4  4  23  2  2  2222 

rOrl  +  rOtlt2r2r3  +  2rOrltlt2r2r3  +  rOtlt2r3 


+  ^rOrltlt2r3  +  VlVs  +  r0rlClr2r3  +  2r0rltlt2r3 


+  r0rlClr2r3  +  2r0tlt2r2r3  +  tlr2r3 


+  6(t  -  Hi) 


55  34  2223  2  23 

r0ri  +  rotit2r2r3  +  2r0tlt2r2r3  +  2r0ritlt2r2r3 


+  ^rQrit1c2r2r3  +  A'AAA  +  ‘vl'lVs  + 


,3  3  2  _  3  2  2  2  2  2  2  3  ,  .  2  2  3 

r0rltlr2r3  +  3rOrltlt2r3  r0ri  lr2r3  +  4r0rlt lt2r2r 3 


.  2  3  2  3  „  3  4  2  24jt  45 

+  r0tlt2r3  +  r0rltlr2r3  +  3r0tlt2r2r3  +  Clr2r3 


+ 


(20) 


Wave-Mode  Coordinate  Solution 


Using  the  wave-mode  coordinate  method  described  in  [2],  the  following 


solution  may  be  obtained  for  F  (t) : 


Fi(t)  (1  -  'o’ 


where 


<5(0  E  (r^)  X((2n+l)x) 

n=0 


CO  oo 

+  E  (r0t2r3tl)n  (  E  (V^'PCn+l,  m)X(2nx)A((2m+l)x)) 


(  E  (r2r3)mp(n,m)X(2nT)X(2mT)) 
'  m=0 


00  00 

+  E  r3t1(rQt2r3t1)n  /  E  (r0r1)mP(n+l,m)A(2nx)X((2-M-l)x)) 

n=0  '  m=0 


OO 

•/  E  (r2r3)mP(n+l,m)X(2nT)X(2(m+l)T)j 
vm=0 


r,/  x  (n  +  in  -  1)  ! 

p(n,m)  (n  _  i)j  m» 


and  X(x)  is  a  time-shift  operator  defined  by 


f (t)X(x)  =  f(t  -  x) 


The  quantities  r^,  r2>  t^  and  t2  are  defined  by  eqns.  (14)  through  (17), 
and  Tq  and  r3  are  again  the  reflection  coefficients  at  the  left-hand  and 


right-hand  boundaries  of  the  lattice,  respectively. 


The  derivation  of  eqn.  (21)  is  discussed  in  detail  in  [2],  Eqn.  (21) 


may  be  obtained  by  setting  x^  =  T2  =  x  (where  x^  and  x2  are  the 
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characteristic  transit  times  defined  previously)  in  eqn.  (103)  of  [2]. 
Writing  out  the  first  few  terms  of  eqn.  (21)  gives 


Fl(t)  ^0  ( R1  +  R2  )  (1  r0> 

•  5 ( t )  A (t)  +  y^Ot)  +  (rQr1)2A(5T)  +...1 

r  2  ~ i 

+  r3t15(t)  A(x)  +  rQr^AOx)  +  (r^.^)  A(5t)  +.  . . 

•  A (2t)  +  r2r^A (4x)  +  (r2r3) 2A (6t)  +. . . 

+  (rQt2r3t1)6(t)  A(3t)  +  2rQr1A(5x)  +  3^^)  2A  (7x)  +. . . 

r  2  _i 

•  A (2t)  +  r2r3A(4x)  +  (r^)  A(6x)  +. . .  J 

+  r3tl(r0t2r3tl)(S(t)  f^'(3T)  +  2r0r1^(5T>  +  3(rQr1)2A(7x)  +.  .  . 

•  J"a(4x)  +  2r2r3A(6x)  +  3(r2r3)2A  (8x)  +. .  . 

+  (rQt2r3t1)26(t)  |a(5x)  +  Sr^A^x)  +  6(rori)  2A  (9x)  +...  j 


Grouping  the  terms  in  eqn.  (24)  according  to  their  time  delays  gives 


Fl(t)  (  R,  +  R  J  (1  "  r0) 


S(t)A(x) [1]  +  6(t)X(3x) !r  r  +  r  1 1 

L°  1  3  !j 


+  <5(t)A(5x) 


(r0rl)  +  r3tir0ri  +  r3tlr2r3  +  r0t2r3tl 


+  5(t)A(7x) 


(r0rL)3  +  r3tl(r2r3)2  +  r^r^r^ 


+  r3t1(r0r1)  +  (r0t2r3tl)r2r3 


+  2(r0t2r3t1)rQr1  +  r3t1(r0t2r3tl^ 


+  6(t)A(9x) 


fcr^r,)4  +  mtJr.rJ3  +  rXxxJr.r 


L 


0  1  3  1  2  3  3  1  0  1  2  3 


+  r3Cl(r0rl)2r2r3  +  r3Cl(r0rl)3 


+  r0t2r3tl(r2r3)  +  2rOt2r3tlrOrlr2r3 


+  3r0t2r3tl(r0rl)  +  2r3tlrOt2r3tlr2r3 


+  2r3tlr0t2r3tlr0ri  +  (r0t2r3tl) 


J 


+  6(t)A (llx) 


<r0ri)5  +  +  rjt  3 


+  r3t1(rQr1)  +  r0t2r3t  1  (-r2r3 


+  2rot2r3tlrOrl(r2r3) 


+  3r0t2r3tl(r0rl)  r2r3 


+  4r0t2r3Cl(r0rl) 


+  3r3C1(rOt2r3tl^r2r3^ 


+  4r3t1(r0t2r3tl^r0rir2r3 


+  3r3t1(r0t2r3ti) (rorl) 


+  2(rQt2r3t1)  r2r3 


+  3(rQt2r3t1)  rQr1 


+  r3tl(r0t2r3tl) 


Comparison  of  Pulse  Summation  Solution  and  Wave-Mode  Coordinate  Solution 


Since  both  eqn.  (7)  and  eqn.  (21)  are  expressions  for  the  same 
physical  quantity  F^t),  eqns.  (7)  and  (21)  give,  when  expanded,  identical 
results.  Note  that  if  point  1  is  a  free  end,  then  r^  =  1,  and  both  eqns. 
(7)  and  (21)  give  F^(t)  =  0,  as  required  by  the  boundary  condition  at  a 
free  end.  It  can  be  seen  from  eqns.  (20)  and  (25)  that  the  first  few 


r  < 1  *!»!'! 


".\VLnW 


terras  of  eqn.  (7)  are  indeed  identical  with  the  corresponding  terms  of 
eqn.  (21) .  Note  that  in  the  pulse  summation  solution  given  by  eqn.  (7) t 
the  impulses  are  grouped  according  to  the  time  delay,  whereas  in  the  wave¬ 
mode  coordinate  solution  given  by  eqn.  (21)  some  manipulation  is  necessary 
before  the  impulses  can  be  grouped  according  to  the  time  delay. 

Equivalent  Paths 

The  numerical  constants  represented  by  the  quantities  N^(K  ,K2>L^) 

and  N^K^.f^.L^)  in  eqn.  (7)  and  the  quantities  P(n+1,  m)  and  P(n,m)  in 

eqn.  (21)  are  due  to  the  existence  of  equivalent  paths  from  the  input 

location  (the  joint)  to  the  output  location  (point  1)  of  the  lattice  in 

Fig.  1.  The  concept  of  equivalent  paths  in  one-dimensional  structures 

consisting  of  multiple  segments  is  discussed  in  detail  in  [1],  Basically, 

a  path  A  from  input  location  to  output  location  is  equivalent  to  a  path 

B  from  input  location  to  output  location  if  a  pulse  which  travels  along 

path  A  arrives  at  the  output  location  with  the  same  time  delay  and  the 

same  amplitude  as  an  identical  pulse  which  travels  along  path  B.  For 

example,  the  four  equivalent  paths  represented  by  the  underlined  term 
2  2  3 

^rgr^t^tj^r^  in  eqns.  (20)  and  (25)  are  shown  on  x-t  diagrams  in  Fig.  2. 
(The  coordinate  x  is  defined  in  Fig.  1.)  An  impulse  of  initial  amplitude 
arrives,  after  following  any  of  the  four  paths  shown  in  Fig.  2,  at 
point  1  with  a  time  delay  of  lit  and  an  amplitude  of 

R1  \  223 

- — — —f  (1  -  r0)rQr1t1t2r2r3  for  each  of  the  four  waves. 

The  pulse  summation  method  is  a  time  domain  method  which  is  based 
upon  a  systematic  enumeration  of  equivalent  paths  within  a  structure  [1]. 


The  wave-mode  coordinate  method,  on  the  other  hand,  is  a  frequency  domain 


method  which  gives  no  explicit  consideration  to  the  existence  of  equivalent 
paths.  The  numerical  constants  which  appear  in  the  wave-mode  coordinate 
solution,  and  which  in  fact  account  for  the  existence  of  equivalent  paths, 
appear  naturally  in  the  wave-mode  coordinate  method  as  a  part  of  the  process 
of  Fourier  inversion  [2]. 


CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  investigation,  nondispersive  pulse  propagation  in  a  one¬ 
dimensional  lattice  structure  consisting  of  two  segments  is  analyzed 
using  the  pulse  summation  method  and  the  wave-mode  coordinate  method. 

It  is  shown  that  the  two  methods  give  identical  results,  and  that  both 
methods  account  for  the  existence  of  equivalent  paths  in  the  lattice 
structure. 

Both  the  pulse  summation  method  and  the  wave-mode  coordinate  method 
may  be  extended  to  nondispersive  pulse  propagation  in  one-dimensional 
lattice  structures  consisting  of  an  arbitrary  number  of  segments.  Such 
an  extension  of  the  pulse  summation  method  is  considered  in  [1], 

Using  the  general  procedures  described  in  [2],  the  extension  of  the 
wave-mode  coordinate  method  to  the  analysis  of  nondispersive  pulse  propa¬ 
gation  in  two  and  three-dimensional  lattice  structures  presents  no  major 
conceptual  difficulties.  The  extension  of  the  pulse  summation  method  to 
two  and  three-dimensional  structures  seems  feasible,  but  has  not  yet  been 
accomplished.  The  exploration  of  equivalent  paths  in  two  and  three- 
dimensional  lattice  structures  may  prove  to  be  very  interesting. 

The  problem  of  dispersive  pulse  propagation  in  lattice  structures 
may  also  be  analyzed  by  the  wave-mode  coordinate  method  described  in  [2] 
with  no  major  conceptual  difficulties.  Computational  difficulties  are 
expected,  however,  for  complex  structures.  The  extension  of  the  pulse 
summation  method  to  dispersive  pulse  propagation  does  appear  to  present 
major  conceptual  difficulties.  In  the  dispersive  problem,  pulses  do  not 
maintain  their  shape  as  they  propagate,  and  it  is  not  clear  how  to  include 
the  effects  of  dispersion  into  the  pulse  summation  method. 
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